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SUMMARY 


A computer program has been developed to prepare wind-tunnel -genera ted 
airfoil data for input into helicopter performance prediction programs. The 
program provides for numerically cross-plotting the data, plotting the data, 
amd then tabulating and punching the tabulated result into computer cards for 
use in the government’s rotorcraft flight simulation model. ' 


INTRODUCTION 

During the process of design and evaluation of rotors and rol-»r airfoils, 
it frequently becomes necessary to use experimentally derived airfoil charac- 
teristic data to predict helicopter performance. Performance predictions are 
usually accomplished through the use of helicopter performance programs such 
as the rotorcraft flight simulation model to calculate forward flight perfor- 
mance and the prescribed wake hover performance program to compute hover per- 
formance. (See refs. 1 and 2.) Generally, helicopter performance programs 
require airfoil characteristic data in a tabular format with angle of attack 
varying with Mach number. These programs generally require a uniform distribu- 
tion of the coefficients at specified values of angle of attack and Mach number. 
This uniformity in distribution of data is rarely achievable in wind-tunnel 
testing; it is more efficient to test through a range of Mach number and angle 
of attack and then interpret the data at selected Mach numbers and angles of 
attack. This approach can result in a process of manually cross-plotting the 
data to prepare them for use in the performance programs. 

The purpose of the current effort is to facilitate the preparation of the 
airfoil data from wind-tunnel experiments and to provide both plots of the data 
and punch card output for use in the performance programs. This purpose has 
been accomplished through the development of a computer program, entitled 
•'PLTAERO,” which is described in this report and is presented in the appendix. 
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SYMBOLS 

blade-element drag coefficient 
blade-element lift coefficient 
blade-element pitching-moment coefficient 
blade-element normal-force coefficient 
Mach number 
Reynolds number 
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angle of attack, degrees 


Oq angle of attack corrected for lift interference effects, degrees 


ANALYSIS 

The PLTAERO conputer program provides for three functions: (1) numeri- 

cally "cross-plotting” the wind-tunnel data to obtain values of c\, Cj, 
and Cgj at specified values of angle of attack and Mach number, (2) plotting 
the data as a function of both Mach number and angle of attack, and (3) tabu- 
lating and punching the data into computer cards for use in the performance 
programs. In accomplishing these functions, the computer program makes use of 
cubic-spline-under-tension subroutines for both the numerical cross-plotting 
and the data plots. References 3 and give the theoretical background for 
these subroutines; the subroutines are unpublished but were vrritten as part of 
the data-reduction computing package for the Langley V/STOL tunr.-l. 


Numerical Cross-Plotting 

The numerical cross-plotting is accomplished by four subroutines, each 
of v.'hich make use of the previously mentioned cubic-sp.'> ine-under-tensiori tech- 
nique to interpolate and extrapolate the input v.’ind-tunnel data. The process 
requires three steps. The first step is to determine whether the wind-tunnel 
data fall within the minimum and maximum values of angle of attack and Mach 
number specified by PLTAERO and its input. The minimum angle of attack is 
required by the program to be and the maximum angle of attack is a func- 
tion of Mach number and is specified by input. The second step in the numerical 
cross-plotting is to calculate the values of the airfoil characteristic data 
at the specified minimum and maximum angles of attack or Mach numbers. If the 
minimums and maximums fall v/ithin the range of the wind-tunnel data, a simple 
interpolation (using cubic splines under tension) is performed. If the minimu-us 
or the maximums I'aH outside the range of wind-tuiinol data, then the cubic- 
spline-under-tcnsion technique is used to calculate the slopes at the minimum 
and maximum points. After the slopes arc calculated, a linear extrapolation 
of the curves is used to obtain values of the airfoil characteristic data at 
the specified "end" points. The third step in the process is an interpolation 
(using cubic sblines under tension) to obtain the data at the specified angles 
of attack and Mach numbers. 

The four subroutines applied are entitled "BOUND!," "B0UND2," "ADJUST!," 
and "ADJUST2.” Subroutine BOUND! calculates the airfoil characteristic data 
at the minimum and maximum angle of attack for each Mach number. Subroutine 
B0UND2 calculates the airfoil characteristic data at the minimum and maximum 
Mach numbers for each angle of attack. Subroutine ADJUST! interpolates the 
data at specified angles of attack for each Mach number and subroutine ADJUST2 
interpolates the data at specified Mach numbers for each angle of attack, the 
first and second steps of the cross-plotting take place in subroutines BOUND! 
and B0UND2; the third step of the process takes place in subroutines ADJUST! 
and ADJUST2. 
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Plotting of Data 


The plotting of the coefficients takes place in two subroutines (PI and 
P2). Plots are presented in two formats. The first format (produced by sub- 
routine P1) presents the curves for cj, c^j, and Cjj, as a function of Mach 
number for constant angle of attack as derived by B0UND1 and ADJUST1 . Fig- 
ures 1 to are sample figures to illustrate the plotting formats of the com- 
puter program. They are all for a Wortmann FX69-H-O98 airfoil with the data 
taken in the Langley 6- by 28-inch transonic tunnel. All the angles of attack 
are presented in one figure, as indicated by the example shown in figure 1. 

The second plotting format (produced by subroutine P2) presents the curves 
for c2, C(j, and Cj. as a function of angle of attack for constant Mach num- 
ber. A separate curve for each Mach number is presented for the blade-element 
lift, drag, and pitching-moment coefficients as shown by the examples in fig- 
ures 2, 3f and Both subroutines make use of cubic splines under tension to 
fair the data. The data presented in figures 1, 2, 3, and are based on the 
wind-turnel data of reference 5. 


Tabulating and Punching the Data 

After the data have been plotted, they are presented in the format illus- 
trated in table I. The information from table I is then "faired” into the 
internal table for an NACA 0012 airfoil of the rotorcraft flight simulation 
model (ref. 1). The table for the 0012 airfoil is included internal to PLTAERO. 
This inclusion is necessary because wind-tunnel tests do not generally 'include 
measurements at angles cf attack from -180° to 180° (i.e., the range of values 
that might be experienced by a helicopter rotor blade as it travels around the 
hub) . 


The fairing is accomplished in the following manner: (1) Because values 

of Cl , c<j, and Cj, at angles of attack less than -4° generally have little 
impact on performance, the coefficients in the reference tables at angles of 
attack less than -4° are left unchanged; (2) the values of C] , c^j, and 
for angles ranging from -4° to otmax maximum specified angle of attack 

at each specified Mach number) arc replaced with corresponding values from 
table I; (3) values of 02, Cj, and Cn, for angles ranging from to 21° 

are interpolated using TSPLII’E; and (4) values of 02, c^j, and 0^ in the ref- 
erence tables for angles of attack ra;ifc.ing from 21° to 180° are left unchanged. 
The resulting table’* may be printed and punched in a format suitable for use 
with the rotorcraft flight simulation program. It is important to note that 
the plots generated by the PLTAERO program can be used to verify that the 
appropriate values have been inserted in the table. Any adjustments desired 
in the tabular data can be easily rectified by comparison with the plots. 


RESULTS 

Figure 5 is taken from reference 5 and shows the variation of C|^ with 
angle of attack corrected for lift interference effects and c<j as a function 
of Cfl. The blade-element lift coefficients used in the present report were 
calculated by using these normal-force coefficients with the corresponding drag 



coefficients of reference 5 to perform an axis rotation to define the lift 
coefficient. Comparison between these lift coefficients and the results of 
PLTAERO is presented in figure 6 as an example of program output. The card 
input guide is presented in table II. 


CONCLUDING REMARKS 

A computer program has been developed to facilitate the preparation of 
airfoil data for input into helicopter performance prediction programs. The 
program provides for numerically cross-plotting the data, plotting the data, 
and then tabulating and punching the tabulated results into computer cards for 
use in the government's rotorcraft flight simulation model. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 17, 1977 
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APPEIIDIX 


COMPUTER PROGRAM TO PREPARE AIRFOIL DATA FOR USE IN 
HELICOPTER PERFORMANCE PROGRAMS 

The computer program PLTAERO is written in the Langley FORTRAN Extended 
Version or FTri ^1, and has been used on the Control Data series 6000 and 
CYBER 175 computer systems under the NOS system. The subroutines and their 
uses are presented in table III. Figure 7 is a diagram of the program struc- 
ture. Basic plotting subroutines, such as those for drawing axes and annota- 
tion, are supplied from the Langley graphics output system; similar routines 
are assumed to be available to the general programing community. PLTAERO 
requires 77400 octal storage locations and takes about 4 seconds to compile. 
Each case takes about 20 seconds to execute on the CYBER 175 computer. 


Subroutine TSPLINE 

Subroutine TSPLINE performs two functions in PLTAERO; it interpolates 
curves at specified values and it calculates the fir.st derivatives of curves 
at specified values. Inputs to TSPLINE include the x- and y-coordinates of 
the curve (x must be strictly increasing), the number of input points, a 
tension parameter (set to 10 in PLTAERO), a computing option, the number of 
points to be interpolated, and the x values at which interpolated values are 
desired. Outputs from TSPLINE include the interpolated values along with their 
first and second derivatives and the area under the curve defined by the input 
X- and y-coordinates. 


Computer Program 

The computer program PLTAERO used to prepare airfoil data for use in heli- 
copter performance programs is presented as follows. 
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TABLE II.- CARD IKPUT GUIDE 


Card 

Variable 

Format 

Description 

I 

1 

NAMAX 

15 

Maximum number of angles of attack 



NMMAX 

15 

Maximum number of Mach numbers 

1 

2 

KA, KB, KC 

315 

Control plotting of Cj, and Cn, respec- 

tively, in PI subroutine; e.g., if KA = 0, 
C| plots arc suppressed, etc. 

1 

! 

3 

LA, LB, LC 

315 

Control plotting of Cp C(j, and c^, respec- 
tively, in P2 subroutine; e.g., if LA = 0, 
Cj plots are suppressed, etc. 

1 

j 

I 

i 

u 

KA, MB, MC 

315 

Control ilculation and punching of tables for 
performance program, e.g.; if MA = 0, 

C| tables are suppressed, etc. 


5 

PMN(I) 

6F10.0 

Desired Mach numbers (extend lowest Mach number 
data back through 0.3, 0.2, to MsO, 

NMMAX values required) 

j 

6 

XNMK(I) 

GF10.0 

Number of desired Mach numbers at each desired 
angle of attack (NAMAX values required) 

1 

e 

f 

7 

PAA(I) 

8F10.0 

Desired angles of attack, using intervals 

starting at (NMMAX values required) 


8 

XWA(I) 

8F10.0 

Number of desired angles of attack at each 
desired Mach number (NM?^AX values required) 

1 i 

i 

! 

9 

m*Ki) 

8F10.0 

Measured Mach numbers (extend lowest Mach num- 
ber data back through 0.3» 0.2, to M r 0, 
NMMAX values) 

1 

\ 

\ 

V 

1 

. i 

10 

XNAT(I) 

8F10.0 

Number of measured angles of attack at each 
measured Mach number (NMMAX values required) 

1 

11 

AA(I,J) 

8F10.0 

Measured angles of attack at Ith Mach number 
(XNAT(I) values) 

1 

1 

12 

CL(I.J) 

8F10.0 

Measured lift coefficients at Ith Mach number 
(XNAT(I) values) 


13 

CD(I,J) 

8F10.0 

Measured drag coefficients at Ith Mach number 
(XNAT(I) values) 


14 

CH(I,J) 

8F10.0 

Measured pitching-moment coefficients at Ith 
Mach number (XNAT(I) values) 

; 

(a) 





15 

Title 

3A10 

Heading for plots (first line) 


16 

Title 

3A10 

Heading for plots (second line) 



*Cards 11, 12, 13, and 14 should be repeated NKMAX times. 
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TABLE III.- SUBROUTINES USED IN PLTAERO 


Subroutine 

Use 

DECKPL 

Input/ output and control 

BOUNDl 

Numerical cross-plotting 

ADJUST 1 

Numerical cross-plotting 

PI 

Plot control 

BOUND2 

Numerical cross-plotting 

ADJUST2 

Numerical cross-plotting 

P2 

Plot control 

C81 

Punch control 

C82 

Numerical fairing of performance tables 

STMBS 

Draws standards, NASA symbol for cubic spline 

CURPLT 

Plots and fairs a cubic spline under tension 
through a set of points 

CUBSPL 

\ 

Calculates cubic spline for CURPLT 

FUNC (function) 

Auxiliary to CUBSPL 

ATAKH (function) 

Auxiliary to CUBSPL 

TSPLINE 

Tension and cubio-spllne interpolating 
subroutine 
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Figure 2.- Blade-element lift coefficient plotted against 
angle of attack at M s 0.30 for FX69-H-098 airfoil, 
generated by subroutine P2. 


1)6 




X 10** 



Figure Blade-element pitching-moment coefficient 
plotted against angle of attack at M s O .30 for 
FX69-H-098 airfoil, generated by subroutine P2, 
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gure 7.- Program structure. 













